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Abstract 

We establish a necessary density criterion for the identifiability of time- 
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condition is based on the density criterion for Gabor frames and Riesz 
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1 Introduction 

The goal in operator identification is to recover an incompletely known operator 
from its action on a single input signal [201 [27]. In mathematical terms: for a 
normed linear space of operators Z mapping a set X into a normed linear space 
Y, we seek g & X such that the induced evaluation map 

^g-.Z^Y, H^Hg, 


1 


2 


is bounded and boundedly invertible on its range. In short, we require that for 
some g € X there exist positive constants A and B with 

A\\H\\z<\\Hg\\Y <B\\H\\z, H G Z. 


This and similar inequalities will be represented by 


||i/5l|y ^ HgZ, 


from now on. 

Identification of operators is important in numerous applications. For ex¬ 
ample, in mobile radio communications it is desirable to identify an unknown 
channel operator prior to information transmission. In radar applications, infor¬ 
mation on a target is gained through analyzing its response to a known sounding 
signal. 

In this paper we focus on Hilbert-Schmidt operators on the space of square 
integrable functions on M, L^(M). Hilbert-Schmidt operators on L^(M) are for¬ 
mally given by 


Hf{x) = jj f{x -t)dvdt 


gnit, v) 7ri(t, v)f{x) dv dt 


with gn G The unitary time-frequency shift vrrf(A) by A = {t,v) G 

is defined by 

7rrf(A)/(x)=TiM,/(x) = e2™("-*)/(x-t), f G (1) 


The space of Hilbert-Schmidt operators HS inherits the Hilbert space structure 
from by setting {H,K)hs = {'nH,gK)L'^ and \\H\\hs = WvhWl'^ [IHlIIl]. 

For A = {s,uj;z,y) G and Hq Hilbert-Schmidt with spreading function go 
we define the operator H\ by 

= Vr2(A)77o = T(^s,oj)M{z,y)V0- (2) 

The central goal of this paper is to establish a density criterion on a not neces¬ 
sarily full-rank lattice A for the identifiability of the closed linear span of 

(//o,A) = {^A}AeA, (3) 

that is, a necessary density condition on A for the existence of g with 


\\Hg\\L2(R) - \\H\\hs, H G span(i?o,A)- (4) 

The paper is structured as follows. Section recalls general facts on modu¬ 
lation spaces, on Gabor Riesz bases and frames, and on Hilbert-Schmidt opera¬ 
tors. In Section]^ we discuss the identification problem outlined above in detail 
and state our main result, Theorem 3.6 Theorem |3.6| is proven in Section]^ 


Section provides some examples of identifiable classes of Hilbert-Schmidt op¬ 
erators. Also, the design of identifiers for the here considered operator families 
is discussed in Section [H 
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2 Background 

This section reviews some basic properties of Gabor Riesz sequences and frames 
in the Hilbert space of square integrable functions and in the space of 

Hilbert-Schmidt operators on 

A countable family of vectors {^aIasA in a Hilbert space H is called a Riesz 
sequence if 

II{ca}||£2(a) II ^ cxgxWn, {ca} e ^^(A), 

AeA 

where f^(A) denotes the space of square summable sequences indexed by A0 
If only ||{ca}||£2(a) > A\\J2\&ACxg\\\n, {ca} e ^^(A), for some positive A, 
then we refer to {fl'AlAeA as Bessel sequence. A Riesz basis is a Riesz sequence 
that spans R. 

A countable family {^aIasA is a frame for T-l if 

ll/llL2(Rd) X ||{(/,5A)}AeA||£2(A)) f (5) 

The set A = C with M being a (not necessarily full rank) real 

2d X 2d matrix is called lattice. A Gabor system {g, A) in is the set of all 

time-frequency shifts 0 of the window function g by elements A = {x,io) G A, 
that is, 

{g, A) = {gx = TTdWg : A G A}. 

The set {g,A) is called Gabor Riesz sequence if it is a Riesz sequence in 
and a Gabor frame if it is a frame for 

Below, we shall utilize modulation space theory as developed by Feichtinger 
and Grochenig [3 0 H E!. Let 5(M) denote the space of Schwarz functions 
on M and 5'(M) its dual of so-called tempered distributions. Let u*, s G M, be 
the polynomial weight function Vs{z) = (1 -|- |^;|)^. The weighted mixed-norm 
sequence space ^s(A) contains all sequences {caIasA with the property that 

ll{cA}||ri(A) = X] 

AeA 

The weighted modulation space (M) consists of all those tempered distribu¬ 
tions / G 5'(M) with 

II/IImi = j \{f,Mzh)vs{z)\ dz < oo, 
with 7 (x) = x G M, and where we refer to 

^ 7 /(^) = (/> ^1(^)7), ^;gM^ 


Recall = L^{A, u) where v is the counting measure on A. 


( 7 ) 
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as short-time Fourier transform of / with respect to the window 7 . 

The weighted modulation space M“(M) consists of all / G 5'(M) with norm 

II/IIm- = sup \{f,TTi{z)'y)vs(z)\ < 00 . ( 8 ) 

2 eK 2 

Note that the dual space of Mf is M“. If s = 0, we write simply and M°°. 
For a detailed treatment of the theory of modulation spaces we refer to Chapters 
11 and 12 of [H]. 

Hilbert-Schmidt operators on L? (M) are in one-to-one correspondence to their 
kernel | 6 l page 267], and, similarly, they can be represented by their time-varying 
impulse response hu, their Kohn-Nirenberg symbol an, and their spreading 
function r]}{- In fact, formally, 


Hf{x) 


j fiH{x,y)f{y)dy = j hH{t,x) f{x -t)dt 
II fix-t)dudt = I 


(9) 


The functions kh, hn, (Xh^Vh are related by 

I 'nH{t,n) du = hH{t,x) = KHix,x -t) = I aH{x,^)e^^^^^ d^, 


and 


= ||/ih||l 2(R2) = ||r//f||L2(K2) = ||fTH||L2(R2). 

due to the unitarity of the L^-Fourier transform H which is densely defined by 


For reference, we include the definition of Beurling density of Let Bd{R) 
denote the ball in centered at 0 with radius R and let |A4| denote the car¬ 
dinality of the set M. The lower and upper Beurling densities of A C are 
given by 


D (A) = lim inf inf 

R-^00 

D'''(A) = lim sup sup 
R -^00 


\An{BdiR) + z}\ 
yolBdiR) ’ 
\An{Bd{R) + z}\ 
vol BdiR) 


In case of D~{A) = D~^{A), we say that A has Beurling density D{A) = D~{A) = 
D^{A). Note that the Beurling density of a lattice A equals the inverse of the 
Lebesgue measure of any measurable fundamental domain of A. See [ 2 T] for a 
more general concept of Beurling density. 
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3 Basic Observations and Main Result 


If {Ho,A) in Q is a Riesz sequence in the space of Hilbert-Schmidt operators, 
then 

span (ffp, A) = {y^cxHx : {ca} G ^^(A)} 

AeA 

and 

II ^ CA-flA||j|^5 - ||{ca}||£2(A), {ca}g^^(A). (10) 

AeA 

In this case, identifiability of span {Hq, A) by g is equivalent to establishing 

II - II{ca}||£2(A), {ca}G^^(A), (11) 

AeA 


that is, to showing that {^Ia^IagA is a Riesz sequence in L^(M). Note that (10) 
corresponds to a Riesz sequence condition in the “large space” L^(M^), while (11) 
is a Riesz sequence condition on a similarly structured family of vectors in the 
“smaller space” L^(M). We shall generally assume that the weaker condition, 
{Hq,A) is a Riesz sequence, holds and focus on the question whether the set 
{Hxg}\£A is a Riesz sequence for some g. 


Remark 3.1 (z) We established that (10) and (0 implies (Q. Moreover, if g 
is square integrable, then 


c\Hxg\\L2 < II ^ cxHxWhs {ca} G ^^(A), 

AeA AeA 


and ([^) can be replaced with the condition that {Hq, A) is a Bessel sequence to 
obtain (^. This argument is not always applicable as we shall generally allow g 
to be a tempered distribution and use the fact that some spaces of operators map 
spaces of distributions into L^(]R). For example, operators in so-called operator 
Paley-Wiener space 


OPW^iS) = {H £ HS{L'^{R)) : snppgn Q A} 


map boundedly the modulation space defined in ([^ to L^(]R) in case 

that S is compact [23]. 

(ii) The identifiability of span {Hq,A) neither implies (10) nor ( [IT| ). Indeed, 
in some cases span {Hq, A) is identifiable, span {Hq, A) = span {Hq, 2A), and 
{Hq,2A) and {RAfl'}AG 2 A are Riesz sequences while {Ho,A) and {Ra^IagA are 
not. 


The framework developed here is motivated in part by the following well 
known results from time-frequency analysis. The first result gives a necessary 
condition on a Gabor system to form a Riesz sequence (see m and references 
therein). 
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Theorem 3.2 Let F = 62 )^^ ® rank lattice in If there exists 

g G L^(M) such that {g,T) = {'7ri(7)5(}^gr is a Riesz sequence in L^(M), then the 
Beurling density D{T) = | det ( 63 )l ^ satisfies D{T) < 1. 


Theorem 3.2 is a special case of the herein established framework of operator 
identification. Indeed, set 


Ar 


A 0 0 0\^ /ai 61 h OV 2 

[0 1 1 0 ) V«2 b2 b2 Oj 


and observe that Lemma 4.1 implies that for any Hilbert-Schmidt operator Hq 
we have 

sp^{Ho, Ar) = spir{7ri(7) o iFo} 7 er- 


So g identifies span (Lfo Ar) if and only if { 7 ri( 7 )(Lfofl')} 7 er is a Riesz sequence. 
It follows that Li(r) < 1 is necessary for the identifiability of span {Hq, Ar)0 
The condition D{T) < 1 implies that Ap is not too dense in the two-dimensional 
“tilted” plane M x {y = z} x {0} C M^. 

The second result motivating this paper plays a critical role in the analysis 
of slowly time-varying operators mM and in the recently developed sampling 
theory for operators 12nil23IMj. 


Theorem 3.3 Let Ai = (^^2 be a full-rank lattice and Hq : L^(M) —: 

A(]R), Hof = p ■ {f * r), be a product-convolution Hilbert-Schmidt operator with 
p,r smooth and compactly supported. If there exists a tempered distribution g 
such that {7ri{'y)HoTTi{'y)*is a Riesz sequence in A(M), then D{M.) < 1. 


Similarly to above, Lemma 4.1 implies that setting 




fO 0 1 
Vo 0 0 



/O 0 ai bi\^ 2 
Vo 0 02 bij 


indicates that Theorem 3.3 may be considered a special case of our framework. 

Theorems 3.2 and 3.3 motivate the central question in this paper which we 
paraphrase as follows. 


Question 3.4 Can we define a density D on lattices A = M7? C so that for 
a positive constant C we have D{A) > C implies span {Hq, A) is not identifiable 
whenever (Hq,A) is a Riesz sequence? 


We choose the following “Beurling-type” density for sets of points A lying 
within general two-dimensional subspaces of M^. 

^Our reasoning uses Hg £ L^(R) and not g £ L^(R). Indeed, the flexibility of choosing 
not square integrable g is quite beneflcial. In fact, for some (Ho, A) we must choose g £ 
M°°(R) \ L^(R) to achieve that {IfAyjAgA is a Riesz sequence in L^(R). 
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Definition 3.5 The “two-dimensional” upper and lower Beurling densities {or 
for short 2-Beurling density) 0 / A C are given by 


Op,(A) 



sup 

t/ew 


lim inf 

R -^00 


inf 

z&U 


\An { 134 {R) + z}\ 

7ri?2 


sup lim sup sup 

UgU R—>-oo z(^U 


\An{Bi{R) + z}\ 

ttR? 


( 12 ) 


where U denotes the set of two-dimensional affine subspaces ofW^ and Bi{R) is 
the centered open ball with radius R inMf. 

If = D^~^(A), then A has uniform 2-Beurling density D( 2 ){A) = 

d- 2 )(a). 

Observe that with 

A=(“‘ 

\a2 02 C 2 d2) 


we have 

^{ 2 ){^) =[(ai&2 — <22^1)^ + (01C2 — a2Ci)^ 

+ (fllC^2 ~ 0 , 2 di)‘^ + (61C2 — 62Ci)^ ( 13 ) 

+ {bid2 - b2di)^ + {cid2 - C2di)‘^] 


Hence, for A_m = 
02 ^ 11 = I det ( 


0 0 ai 61 
0 0 02 62 
oi 02 M -1 

bl 62 ) I 


)^Z2 in Theorem 
, and for Ap = ( 


3.3 


1-1 


ai bi bi 

a2 b2 62 0 


we have D(^ 2 ){^m) = |ai ^2 - 
^ in Theorem 




3.2 


we have 


D(2)(Ar) = 2-V2|det(?(^^^)|- 

The main result in this paper provides a necessary condition on A for the 
existence of g so that {Hxg}\^K is a Riesz sequence. Mf and M°° denote 
modulation spaces whose definitions we recalled in Section 


Theorem 3.6 Let A = Ml? C and Hq be an operator with rjHo £ Af^^(M), 
s > 2. If {IIq,A) is Riesz and span {Hq, A) is identifiable by g € M“(]R) then 
71(2) (A) < \f2. 


If the hrst or the last row of M is zero, then the bound on 11(2) (A) can be 
improved. Identifiability then implies 11(2) (A) < 1. This is the case in Theorem 


is not clear whether the constant \/2 in Theorem |3.6| can be improved in the 
general case. 

Identifiability of span {Hq, A) depends on both, Hq and the lattice A, so it is 
not surprising that there exists no c > 0 with I1(^^(A) < c implies span {Hq, A) 
is identifiable. In fact, it is easy to construct for any e > 0 a Riesz sequence 


and [13] I. 


{Ho,A) with Z1(2)(A) < e but span(i7o5A) is not identifiable (see Example 


5.4 


3.3 


and Theorem 


3.2 


respectively (in Theorem 3.2 even 11(2) (A) < 2 ^?). It 
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We would like to emphasize that the operator outputs considered herein are 
in the “small” space L^(M) while the kernels and spreading functions of the 
operators are in the “larger” space of functions. This dimension mis¬ 

match implies that a single evaluation map H ^ Hg cannot identify the space 
of Hilbert-Schmidt operators as a whole, just as a degree of freedom counting 
argument shows that the space of complex n x n matrices requires the use of n 
identifiers for identification. 


4 


Proof of Theorem 


3.6 


The following lemmas are used in the proof of Theorem 3.6 


Lemma 4.1 Let go denote the spreading function of Hq G Then the 

operator TaM^T-cHoTcMd has spreading function 


VTaMiT-^HoT^Md, = Tafi+dMb^cVO, b,c,d € M. 


Hence, 

hH= cxVHx = Y CAvr2(A)r/o. 

converges in -norm if and only if 

H = Y^ cxTsM^T_yHoTyM^_^ 

X={s,uj,^,y)^A 

converges in HS{M). 

Proof. A change of variables t = t — a, v = u — b — d and the relation 
TxAdu] = M^Tx implies 


Ta,b+dMb,crio{t, i^)fix - dvdt = 

^ 2 ni{b{t-a)+c{u-b-d))^^^^ - a, i/ - 6 - d)f{x - dndt 

^2ni{bt+cu)^^^^^ i/)/(x - t - dn dt 

i^)Tt+aMj,+6+rf/(x) dudt 


= TaAIbT-c JJ go{t,v)TtMyTcMdf{x)dvdt 
= TaMbT_,HoT,Mdf{x). 
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Hence, in particular, 

{Hf ){x)= JJ riH{t,v)f{x dudt 

= ^ cxffTs^uiM^^yr]o{t,u)f{x-t)e^^"''^^~^'^dudt 

\={s,uj,^,y)&A 

= ^ cxTsM^T_yHoTyM^_^f{x). 

X=(s,oj,^,y)^A 


□ 


Lemma 4.2 Let p,q G C'“(M'^) and let P be the product-convolution operator 
with spreading funetion gp = p® q. Then there exist functions V'l, V ’2 £ 5(M'^) 
such that 


\Pf{x)\ < ||/||m=o(m.)|V’i(^)I, I WMI < ||/||m»(m^)I^ 2 M|, / G 


Lemma 4.2 is a straightforward generalization of Lemma 3.4 in [20] and its 
proof is omitted. Lemma |4.2| with d = 1 is used to prove the following result. 


Lemma 4.3 For H with spreading function po G s > 2, there exist 

pi{t), p 2 {t) = 0{t~^) with 


\Hf{x)\ < ||/||m-(r)V^i(^), \PHfiio)\ < WfWM^mM^), f e M“(M). 

The decay estimates show that Hg ^ L\R) HH (2.52)]. 

Proof. By choosing parameters a,b,c,d with ac < l,bd < 1 and functions 
p,q€ C“(M) with [-§,§] C suppp C [- 5 , 3 ] ^ suppg C [-^,^] we 

obtain a Gabor frame {pp, aZ x 6Z x cZ x dZ) for where pp = p®q (see 

(23 pages 22-23] or |2^). Because pp G 5(M^) C M^^(M^), the Gabor system 
{pp, aZ X bZ X cZ X dZ) is a universal Banach frame according to the definition 
of Grochenig for all modulation spaces (m 2) [a Section 13.6]. The main 
result from [T6| states that the canonical dual window satisfies pp G Mj^(M^). 
Let us denote by the weighted mixed-norm sequence space containing all 

sequences such that the norm 

Il{cj||ri(z 4 ) = \czVs{z)\ < 00, 


where Vg is defined by (|^. 

Then the result from m Corollary 12.2.6] shows that the series expansion 
/= Y/ {f^'Pak,bl^cm,dnhP)Tak,blP^cm,dn'>lP^ f^Mgi^'^) (l^) 
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is convergent in the M^-norm and 

II/IImi(R) - \\{{f,Tak,bl^cm,dnfip)}\\el{Z‘i), f ^ 


(15) 


Furthermore, as the coefficients have decay stronger than the convergence 

of the series holds in L^{M?) and in 

Since the operator H has spreading function 'qn £ (14) and (15) 

imply 

'HH — ^ ^ (^k,l,ra,n^ak,bl^cm,dn'nP 

k^l,m,n£7j 

for some {ck,i,m,n} £ It is legitimate to use as identifier distributions 

g G because of the inclusion 


5(M) C M^(R) C M^(M) C M°°{R) C Ml 


cS'i 


This inclusion is a consequence of im Proposition 11.3.4] and [HI Corollary 
12.1.10]. The fact that the constant weight 1 is n^-moderate ( |14[ Lemma 11.1.1]) 
provides the inclusion Mg{R) C M^(M). 

Next, we estimate the decay of in the time and frequency domain. We 
shall use the fact that translation and modulation are isometries on M^(M) and, 
hence, also on M°°(M) and estimate 


iff9Wi = i E TakM- 

cm T—dnPTdn^bl—cm 5 '(^)| 

— ^ ^ ‘Pak—dn\P'Pdn^hl—cm9{sP)\ 

k^l^m^n^'L 


(16) 


where (pi 

\\9\\m- 


4.2 


For the sake of clarity we denote the expression 


G 5(M) is some positive function, such that 
i(x) after Lemma 


on the right-hand side of (lOiby 


^i(®) = 

k^l^m^n^.'L 


I 'Pak—dn4^l{^^ — E ^k,nTak— dn4^1 (^) 1 

k^n^lj 


where Ck^n — 

We claim that $i(x) = 0{x~^) due to {\ck^i^m,n\}k,i,m,nGZ ^ ^ 

Let us make a change of variables |xf = y > 0, which is equivalent to x = y^ ^ 
2 

x > 0, and X = —y^ , x < 0. Then 

sup |x^| • I E ^k,nPak—dn4^l{^')\ I?/ ^ ^ (?/ ^ )| • 

fc,nGZ k,neZ 
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Since is monotone on [0,oo) and due to our choice (pi G 5(M) (that is, (pi 

~ 2 

decays faster than the reciprocal of any polynomial on M), p[y) = pi{y^) also 
decays faster than the reciprocal of any polynomial. Then we can estimate 
snVy>o\y^Y.k,n^k,nTak-dnP{y)\ by Setting ^'i(y) = yp{y),^2{y) = y‘^^{y) and 
using the equality y"^ = {y — ak — dn)'^ — {ak — dn)'^ + 2y{ak — dn). Hence, by 
the triangle inequality it follows for 0 that 


y ^ ^ ^k,nTak—dn^{y) 
k^n^’L 


^ I ^ ^ ^k^nTak—dn^2 
I ^ ^ ^k,ni^k dfl^ 

21 ^ ^ dfl^Tak—dn^ 1 


(18) 


Taking the supremum in (18) leads to 


sup 

y>o 


y ^ ^ ^k,n'^ak—dnPiy') 
k,nGl^ 


< sup 
y>o 


^ ^ ^k,n^c 


ak—dn^2 


k,n 


+ sup ^ Ck,n{ak - dnfTak —dn^ 


y>o 


k.n£'L 


(19) 


+ 2sup ^ Ck,n{ak - dn)Tak-dn^i 


y>o 


k,n£Z 


We compute the following upper estimate of the summands on the right-hand 


side in (19) 


sup 

y>o 


X] Ck,nTak-dn<Piy) < ^ || SUp |T2(y) | 

k,n&Z k,n&Z 


+ ^ |cfc,n(aA; - (in)^l sup |(^(y)| (' 2 Q') 

y>0 '' ^ 


k^n^’L 


-F2 ^ |cfc,n(afc - dn)| sup |Ti(y)|. 


k.n£Z 


y>o 


We analyze separately the three summands from (20). Since (/>, Ti, ^'2 belong 
to the Schwarz class, they are bounded and decay faster than the reciprocal 
of any polynomial. Also the fact that {ck,i,m,n}k,i,m,n&z £ ■^s(^^) implies the 
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existance of constants Ci,C 2 > 0 such that 


\ck,n\ ■ \a.k - dn\‘^ < Cl ^ |cfc,„|(l + a|A:| + d|n|)^ 
k^n&, 

< ll{Cfc,Z ,m,n}fe,Z,m,nSz||£l(Z4) ^ CX3, 

y~] |cfc,n| • |aA: - dn| < 6*2 ^ |cfc,„|(l + a|A:| + <i|n|) 

k^n^lA k,n^l^ 

< \\{Ck,l 


Furthermore, 


|cA;,n| < l|{cZc,Z ,m,n} k,l,m,n(lZ 

k,n£Z 


< OO. 


Thus the expression on the right-hand side of (20) is bounded, implying the 
desired decay rate oi Hg for a: > 0. In a similar fashion we prove the decay for 
X < 0. Thus sup^-gK |3z^$i(x)| < C and \Hg{x)\ < ||fl'||M°°(R)V^i(3^) has decay 
0(x“®), s > 2. 

A similar estimate can be done for the decay of the Fourier transform ol Hg 


i^«9MI = i E Ck,l,m^nM—akT—crnM(inJ'PT(inMi,i—c^g(^Lij) \ 
fc,Z,m,nEZ 

^ FI \Ck,l,m,n\ • T-cm\J~PTdnMbl —cm 
k,l,m^n£lA 

< Hull M°°(R) ^ ^ \ck,l,m,n\ ' T—cm4’2{^)i 

k,l,m,n^Z 


( 21 ) 


where 4>2 £ 5(M) is some positive function, such that \PPTdnMf,i-cmgi^)\ < 
llfflljV^oo(R)(/> 2 (az) after Lemma [T^ We denote the expression on the right-hand 
side of (21), by 

*^* 2 ( 3 ^) = ^ ^ \Ck,l,m,n\P—cm4^2{x), 

k,l^m,n 

and prove in a similar fashion that d> 2 (a^) = 0{x~^), s > 2. □ 

Proof of Theorem 3.6 For m, n G Z and A = M{m, n)^, we observe that 


Px P'(ai—di)m+{a2—d2)nP^cim+C2r 


PoPdim+d2nPd^bi —ci )m+(b2—C2)n • 


Set gm,n = HoTdj^rn+d 2 nM(^bi-ci)m+{b 2 -c 2 )ng- Since G Mf{R), s > 2, Lemma 
implies the existence of fiix) = 0{x~‘^), 4 > 2 i^) = 0{u}~^),s > 2 such that 


4.3 


\gm,n{x)\ < </’l(3^)||fi'llM°°(R)) 

\Pgm,ni^)\ < M^)\\g\\M-°(R)- 


( 22 ) 
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To prove the claim of the theorem, we show that the family 


{'^{ai—di)m+{a2—d2)n^cim+C2n9m,n}m,n£'Z^ — ^ 

is a set of Gabor molecules [3l Definition 10]) localized with respect to the lattice 

^ai — di 02 — (1-2^ 


A = 


Cl 


C2 


Z" 


For that purpose it suffices to show the existence of T G W{C,£‘^), such that 
K 5 'm,n,TrM^ 7 )| < T(x,a;) for all (m,n) G G Here -f{t) = 

and W{C,£^) denotes the Wiener amalgam space consisting of all continuous 
functions / on such that the norm 


\\f\\w(C,P) 


Note (22) implies 



< oo. 


\{9 T.M^-f)\ < {\9 m^n I ? Tx ItI) = \9m,n\ *l{x) 

< ll5llAf-(R) • </>l *l{x) , 

I (-T 9m, M.,T^-f)\ < {\g 

m^n \,T^h\) = \9m,n\ *1(1^) 

< llfi'llAf°°(R) • 4>2 * 7(‘^)) 


see [I5l[22]. Hence, by setting 


Kt) = II5||m<-(m) max{(?;>i =t= 7(t), =t= 7(-t), (1)2 * lit), </>2 * 7(-i)}, 

we obtain 

1(5 m^ni Tx Afa;7)l — h(max{|x|, |w|}) — ^(ll(^)^)lloo) 
with \h{t)\ = s > 2. Hence, there exists a constant c such that 


\{ 9 m,n,TxM^j)\ < c - h{\\{x,uj) 


that can be bounded in turn pointwise by a function T(x,a;) G W{C,i‘^). Thus, 
{Hx9}\£A is a set of Gabor molecules. 

Assume that span (Hq, A) is identifiable by 5 G M“(]R). By employing [3l 
Theorem 8 ] and [21 Theorem 3] (the latter result is a Gabor molecule extension 
of Theorem 3.2) we obtain that the Beurling density of A, given by 


D{A) = 


det 


ai — di 02 — d2 
Cl C2 


-1 


must be less than 1 in order for {Hxg}\^\ to be Riesz in L^(M). A computa¬ 
tion shows that this condition follows from D(A) > x/2, yielding the bound in 
Theorem 13.61 □ 
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Remark 4.4 To obtain similarly a necessary density condition for to 

be a frame, we would have to show that |(ai — di)c 2 — (02 — d 2 )ci| = < 1 

follows from D{A) > c for some positive constant c. But this is not reasonable 
to expect, as increasing bi and/or 62 greatly in 

D{A) =[(0162 - 02 ^ 1 )^ + (aiC 2 - a 2 Ci)^ 

+ (oi<^2 ~ a2(ii)^ + (61C2 — 62 Ci)^ 

+ {bid2 - b2di)‘^ + {cid2 - C2di)^]~^^^. 

would decrease D{A). 

5 Examples of identifiable operator classes, design of 
identifiers 

To establish identifiability of span {Hq, A) we seek an identifier g such that any 
choice of coefficients {ca} G £^(A) in H = can be computed from 

Hg. Equivalently, we require that {ca} can be computed from the values of the 
Gabor coefficients = {Hg, 7ri{fi)'y), g G At, where ( 7 , At) is an L^-Gabor 
frame for appropriately chosen window 7 G L^(M) and lattice At C Then, 
our task is to solve the system of linear equations 

VfM = {Hg, = '^cx {Hxg, -n-iigh) = ^ A^;A ca, /X G At 

AeA AeA 

for {ca}. The doubly infinite matrix A has entries A^;A = {Hxg,TTi{g)^). 

If g is such that the map A: i‘^{A) —>• has a bounded left inverse then 

span {Hq,A) is identifiable. 

The design of identifiers g can be carried out on the coefficient level. In fact, 
when ( 7 , At) is an appropriately chosen Gabor frame for L^(M), or, for example, 
an £“-frame for M“(M) [T], then we seek a coefficient sequence {d^} so that 
the bi-infinite matrix with entries 

A^-,x = dj:{HxTTiig)j,7Ti{gh) 

JleM 


is invertible. 

To illustrate the method outlined above, we give an alternative proof of one 
implication in Theorem 3.1 in m- Also see [21] for a comprehensive treatment 
of sampling and identification in operator Paley-Wiener spaces. 

Theorem 5.1 The operator Paley-Wiener spaee OPW‘^{[0,a]x[0, ^]) is identi¬ 
fiable. 
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Proof. By definition the given operator space consists of operators H G 
HSfM.) such that 


r]H G span{Mi^^^^X[o,a]x[o,j^]’ kjGZ}, 


that is, we consider 


A 


/O 0 0 
Vo 0 a 



Z^. 


Since : {x,y) G [0, a] x [0, ^],/c, £ G Z} forms an orthonormal basis 

for the space L^([0, a]x[0, ^]) [l9], the spreading function of any operator H G 
OP1B([0, a] X [0, ^]) has a unique expansion 


VH = '^ Ck/Mk^^f]Q 
k,i&Z 


with r]o{x,ui) = X[o,a](a;)X[o,i](‘^)- 

Set Ck^e,m,n = Ck^i6ofi{m,n). We choose 7 = a~2X[Q^a] observe that the 
Gabor system ( 7 , aZ x is an orthonormal basis for L^(M). Using the formal 
relationship {Hg, f)= {'nH,ygf)L2(R) [HI Lemma 3.2] with g = 5az = 
J2n&z^na G M“(M][^we Compute 

■^p,q',k,l {kik,i^aZi ^^'yaq'y'} 

~ 2 A/fe_p ^ 70 , 14 ^ 27 ). 

The Zak transform Za satisfies the relations 

_i _i 

^X[0,a] ^X[0,a] 

and 

a-^^aX[o,a](a:,a;) = e2™[S]-, 

for which we refer to m Chapter 8 ]. Then 


A 


p,q;k/ 


( I 27r2f^— p)(a;+aq) 2'Ki(u}-\-^\ai —2'Kia\ — \u 

X[o,a]{x + aq)e Xro idw + ^je e 


dx du 


We make the substitutions y = x + aq, z = uj P ^ and note that since the 
integrand is nonzero for aq < x < aq + a, [|] = q, it follows that 


A. 


■p,g;k,e — I / e 
Jo Jo 


“ 2m— — 2malz—2'iTiaq(z—-) 

P a P a J 


dy dz 


^^2m(pq) r e^^^^^^dyx /“ 


= Sp,q{kJ), 


^Note that the inner product is still well-defined since rjH has compact support. 
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where we used that : n G Z} and : m G Z} are orthonormal 

bases for L^[0,a] and L^[0, ^], respectively. 

The matrix A = {Ap^q-k/)p,q;k,iei is the identity, and thns invertible, which 
is what we had to prove. □ 

Recall that Vgf{z) = {f,TTi{z)g), z G is short-time Fourier transform of 
f ^ L? (M) with respect to the window g ^ (M) m Chapter 3]. An additional 

positive identifiability result is the following. 

Proposition 5.2 Let /i, {^aIasA £ L^(M) be such that {^'a} is a Riesz se¬ 
quence for its closed linear span in L^, ||/i||i 2 (-]R) = 1. Then the operator family 
span {H\, gHx = ^hg\\, is identifiable. 


Proof The Riesz sequence property of {sja} implies that {gH^} is a Riesz 
sequence as well since it is the image of a Riesz sequence under the unitary map 
Vh ■ 9x ^ VHx 0. Consider the associated action of the operator Hx on g: 
H\g = (7 a (s') h). Therefore, the entries of the matrix A have the form 


Ag-x = dix{HxTTi{ft)j,gp) = dg{gx,gp){7ri{il)g,h). 






(23) 


Whenever { 7 /^} is chosen biorthogonal to {sa}) (23) becomes 
Ap-x = 5{iJ,-\) Y c^A(^i(Ai)7) ^)) 


which shows that the matrix A is diagonal with non-zero entries for appropriate 
{dp}. Hence, A is invertible. □ 

The following examples address more involved rank-2 lattices, for simplicity 
of calcnlation, we will consider Gaussian kernels only. 

Example 5.3 Let Ho be given by ko{x,oj) = ^ that is, go{t,n) = 

J_ „-7riC2iv„-f (F+i/2) 

1. Let A=(Q|g)(g)^Z^. If a, ft are such that \a{fl + a\/2)\ > \/2, |a,0| > 
a/2, |a| > 1, then the operator family span(Ffo: A) is identifiable. 

2. Let A = (ooq^)^Z)^. Ifct,ft are such that |(a| > 1 and G Q, then 
the operator family span {Hq, A) is identifiable. 


The pairs {a, ft) satisfying the conditions in Example |5.3[ Statement 1, are 
illustrated in Figure Statement 2 of Example |5.3| indicates that identifiability 
may depend on conditions that are not expressible in terms of density or simple 
inequalities. When ^ Q, the Gabor system {go, A) is not a Riesz sequence, 
so the problem is not considered in this paper. 

The following example shows that D{A) being small does not necessarily 
guarantee identifiability of span {Hq, A). 
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Figure 1: The set (a,/3) fulfilling the conditions in Example 5.3, Statement 1, 
is represented by the shaded region. The region encompasses the intersection of 
the convex hull of the parabolas |a(/3 + V^a)] > V2, \af3\ > 2 with the subset 
of the plane |q;| > 1 . 


Example 5.4 Let Hq be given by r]Q € s > 2, and let A = 7?. 

/f |a/3| < 1, then the operator family span {Hq, A) is not identifiable. 

The condition \a/3\ < 1 cannot be expressed in terms of 2-Beurling density of 
the index set A, D(Ai) = - f In fact, for any e > 0, we can find a, fi with 

\aj3\ < 1 such that D{A) < e. For instance, when a = 10^°,/3 = (10^° + 1)“^, 
|a/3| < 1, so the family span [Hq, A) is not identifiable, but D[K) k, 10“^*^ is very 
small. 

Further examples of identification using the approach described in this paper 
are given in m Sections 4.4 - 4.6]. 
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